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Abstract 

We consider the 2D critical Ising model with spatially periodic boundary conditions. It is 
shown that for a suitable reparametrization of the known Boltzmann weights the transfer 
matrix becomes a polynomial in the variable csc(4w), being u the spectral parameter. The 
coefficients of this polynomial are decomposed on the periodic Temperley-Lieb Algebra by 
introducing a lattice version of the Local Integrals of Motion. 

1 Introduction 

In this paper we push our investigation of the critical Ising model a bit farther. Our objec- 
tive is to extend the results of the previous work on Critical Dense Polymers [6] to the Ising 
work, in order to put the previous work on a more solid ground and to show that once more 
the same approach is not only possible, but can in this case be carried out completely. 
The approach goes all the way back to [3], where the integrable structure of conformal field 
theory (CFT) was analyzed from the point of view of the existence of infinitely many higher 
spin conserved quantities in mutual involution. Inspired by this work we seek and succeed 
to obtain the complete tower of commuting conserved quantities. 

The involutive charges of CFT are defined as the maximal Abelian subalgebra of the envelop- 
ing algebra UVir of the Virasoro algebra, the lattice analogue of this structure is obtained 
by replacing the Virasoro algebra with the Temperley-Lieb algebra, as it is well known [7J. 
There are of course different types of Temperley-Lieb (TL) algebras, in particular it makes 
a lot of difference to work with open boundaries or periodic boundaries. In j6] we worked 
with an open boundary TL algebra, and obtained closed form expressions for many of the 
involutive charges, an ansatz closed form for the generic expressions was also provided in 
which some of the coefficients are still unfortunately not known. The heart of the difficulty, 
as pointed out long ago also in [7J lies in the presence of boundary terms which manifest 
themselves through the appearence of boundary tangles. It is therefore natural to expect, 
and this is indeed the case, that such difficulties can be overcome by considering a model 
based on a periodic TL algebra. In this case there are no boundary tangles and the ex- 
pressions for the involutive quantities simplify enough that a closed form for them can be 
found. This is surely the case for Critical Dense Polymers on the Cylinder [§] as well, since 
it shares an extremely large number of analogies with the Ising model among which are 
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almost identical expressions for inversion identities, eigenvalues of the transfer matrix and 
expression for the continuum involutive charges. 



2 Generalities 



It is well known that the Ising model can be characterized in terms of Boltzmann weights, 
which can be derived from an elementary lattice hamiltonian [1] . Naturally such weights can 
be suitably rescaled by multiplying for trigonometric funtions, in particular we shall employ 
here the following parametrization: 
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A transfer matrix can then be defined in different ways and for different boundary conditions, 
in particular if we choose to work on the cylinder we can build the transfer matrix in a 
boundary condition dependent way evolving into periodic imaginary time, or on the contrary 
we can choose to build a boundary indipendent transfer matrix with spatially periodic 
boundaries which evolves in non compact time. We choose here to explore the second 
possibility, the first one having been already exhaustively discussed in [Ij and also elsewhere. 
Therefore we introduce the periodic transfer matrix: 



t n=l 
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with 



T, 



n+L 



and similar conditions for a, p. The transfer matrix has thus size 2 L x 2 L . 
The partition function is then defined as: 



TrT 



(2.4) 



(2.5) 



(2.6) 



on general grounds the partition function is expected to behave asymptotically as a linear 
combination of CFT characters, whose modular parameter is related to the exponential of 
some trigonometric function of the spectral parameter. We will describe this behaviour of 
the partition function in a later section. 

It is also well known that the Ising model and in generally Spin-chains, Potts models, RSOS 
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models, Vertex models [7] are directly related to representations of the Temperley-Lieb (TL) 
algebra. We will extensively investigate this connection for this special case, thus for our 
needs it is necessary to introduce the periodic TL algebra, whose generators are defined by 
the following matrix elements: 

L 

(e 2 i)a, P = V25- ai!ai+1 Y[ 6 ak , Pk (2.7) 
k=i 



{e 2 i-i)a, P = -j= n 5* k , Pk (2.1 



which satisfy the relations: 



,2 



&• = V2ei (2.9) 
eie i±l ei = (2.10) 
e^ej = ejei , \i — j\ > 2 (2-11) 

3 Inversion Identity and Eigenvalues 

In this setion we introduce a quadratic functional equation satisfied by the transfer matrix 
T, the techniques here employed naturally follow the work pQ, needless to say. We call it 
Inversion identity, even though strictly speaking its right hand side is not proportional to 
the identity, nonetheless this is the historical name introduced by Baxter in his book [10] in 
a very similar case and therefore we shall use it. 

We also drop for the rest of the paper the use of trigonometric functions of the spectral 
parameter u in favour of the extremely more convenient variable: 

x{u) = sin(4w) (3.1) 

notice that 

x (u+ ^) = -x(u) (3.2) 
One has then that the transfer matrix satisfies the following functional equation: 

T(x)T(-x) = (-l) L 2 2L+1 (T 2i (l/x)l + R) (3.3) 

where T n (z) are the Chebyshev polynomials of the 1st kind, defined by: 

T n (cos#) = cos(n9) (3.4) 

and R is an involution 

R 2 = 1 (3.5) 
which also commutes with the transfer matrix T: 

[T(x),R] = (3.6) 
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it is thus possible to diagonalize the involution and the transfer matrix simultaneously, and, 
calling r = ±1 the eigenvalue of the involution, and A(x) those of the transfer matrix, one 
has: 

A(x)A(-x) = (-l) L 2 2L+1 (T 2L (l/x)+r) (3.7) 
which can be factorized as: 



A(x)A(-x) =(-l)*2" | | ( ^ + sin ( 2fc ~^ L 2Mod(L ' 2) 7r 



2L / 
k=l \ 



(3. 



which sharing out the zeroes leads to the following factorized form of the eigenvalues: 

A(s, /*, r, L) = ^ H (l - xx fe z sin ~ ^^ tt) j (3.9) 

where: 

xx = (/xi, . .. (3.10) 

A*2 = l (3- 11 ) 
given /i we can also introduce the associated Fermionic Partition V = (hi, . . . , k p ) (p < L), 
defined by: 

k G V, if fi k = 1 (3.12) 

for example /i = (1, —1, —1, —1, 1, 1, —1, 1) is associated to the partition V = (1, 5, 6, 8). 
Now, by calling Aq the eigenvalue corresponding to /x = (1,1, 1,1,1,.. .),and introducing the 
auxiliary functions: 

1 + xsin ((2k - <5 r ,i_ 2 Mod(L,2))^ ) 
Pl(x) = ) (- (3.13) 

1 - ISm [(2k - 5 r> i_2Mod(L,2))^J 

one has that the generic eigenvalue is decomposed as: 

A r (x) = A r (x)l[pl(x) (3.14) 

we finally remark that since the multiplicity of the solutions of the inversion identity is 2 L+1 
we have to impose the following selection rule to give account of the actual dimensionality 
of the transfer matrix, which is only 2 L : 



L 

r = 

k=l 



I/x fc = (-l) L ^ (3.15) 



We also introduce the following abbreviation, which will be used whenever its meaning is 
not ambiguous: 

£r,l-2Mod(L,2) = $ (3.16) 
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notice that, for even L 

5 = 8 rA (3.17) 

whereas for odd L 

5 = 5 r ,_i (3-18) 

4 Temperley Lieb Decomposition 

In this section we are going to decompose the transfer matrix T on the TL algebra. This 
task has carried out by analyzing the decomposition with Mathematica up to L = 9, the 
general behaviour has been then recognized to hold for generic L. 
We will decompose the tranfer matrix according to the following expansion: 

T W = ^rE D ^ ( 4 - 1 ) 

n=0 

we then introduce: 

2L 

[e n+ 2, • • • [e n+ 2fc-3, e n+2fc-2 ] • • •]]] (4-2) 

n=l 

A 2 = -2L1 (4.3) 

A 4 = -36L1 (4.4) 

A 6 = -2400L1 (4.5) 

A 8 = -352800L1 (4.6) 

Aio = -91445760L1 (4.7) 

Ai = Hi - 2L1 (4.8) 

A 3 = H 3 + 6Ai (4.9) 

A 5 = 6H 5 + 60A 3 - 120Ai (4.10) 

A 7 = 90H 7 + 210A 5 - 5040A 3 + 5040Ai (4.11) 

A 9 = 2520H9 + 504A 7 - 45360A 5 + 604800A 3 - 362880A! (4.12) 
A n = 113400H n + 990A 9 - 221760A 7 + 11642400A 5 - 99792000A 3 + 39916800A! (4.13) 
and in general: 

_ (2fc-2)l ^ (_i)m+i / 2fc - m - 2 \ (2fc-l)l 

A2k ~ l ~ ^^ Ii2fc - 1 + ~^r~ \ m-l )(2k-2m-2)l A2m - 1 (4U) 



the first D n then read: 



D = 1 (4.15) 
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Di = ^Ax (4.16) 

D 2 = i(A? + A 2 ) (4.17) 

D 3 = ^(A 3 + A? + 3A 2 Ai) (4.18) 

D 4 = -J-(4A 1 A 3 + A 4 , + 6A 2 A? + 3A* + A 4 ) (4.19) 

D 5 = ^r?{Al + 10A?A 2 + 15AiA* + 10A?A 3 + 10A 2 A 3 + 5AiA 4 + A 5 ) (4.20) 

D L = ^((-1) L 1 + R) (4.21) 

Dn = ^! 5 "( Al '---' A ^ ^ 

where the B n (x\, . . . ,x n ) are the Bell polynomials, which are generated by the following 
expansion: 

«P ( £ ^f) = 1 + t ^p^t" (4.23) 

n=l n=l 

We call the A 2n _! the Lattice Local Integrals of Motion (IOM), since they play a completely 
analogous role as the I 2n _i of [3], as we shall realize even better in the next section. 
The A 2n _! form a maximal Abelian subalgebra (of hermitean operators) of the enveloping 
algebra of the TL algebra, in the same way as the CFT Local IOM I 2n _i span the maximal 
abelian subalgebra (of hermitean operators) of the enveloping algebra of the Virasoro alge- 
bra. As remarked in [7] the 2 enveloping algebras should become isomorphic in a suitably 
defined L — > oo limit. 



5 Local Integrals of Motion 

In this section we are going to derive the exact expressions for the eigenvalues of the lattice 
Integrals of Motion and after this we will discuss how these eigenvalues are related to the 
eigenvalues of the conntinuum Integrals of Motion which have been already obtained by 
means of Themodynamic Bethe Ansatz (TBA) in [5], following lines advocated in [2] • In 
this case however no TBA will be employed, rather the eigenvalues will be recognized in an 
expansion that has been already carried out for Critical Dense Polymers in j2] by means of 
all order Euler-Maclaurin expansion of the factorized form of the eigenvalues. In this case 
however we will not employ this method, but rather direct Taylor expansion. 
We start by reminding that 

(2k -5) , . 

tk = 2L 71 (5 - 1} 
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one then sees that the eigenvalues of the transfer matrix can be expanded in the following 
series: 



x 2n 



log A r (x) = -L log x + 2L log 2 - ^ ^T^ry (2n - 1)!2 2 " ^ sin 2 " t fc + 

n=l ^ k=l 



oo 2n— 1 

+ E rarrw t 2 " - 2 ' !22 ""' ( 2 E ( * - E ■">*-' ( * 

n=l fee?? fc=l 



(5.2) 



we then recognize the eigenvalues of the operators A n : 

L 

A 2n = -(2n-l)\2 2n Y,^ 2n t k (5.3) 

k=l 

and 

L 

A 2n . 1 = (2n - 2)\2 2n ~ 1 (2j2^ 2n ~ 1 tk - ^ sin2n_1 **) ( 5 - 4 ) 

keV k=l 

these values of the eigenvalues are generally correct, what happens is that for L > n — 
1 the numerical diagonalization of the Temperley Lieb decomposition of A 2n _i given in 
the previous section agrees perfectly with the above expressions, below this threshold the 
Temperley Lieb decompositions of the previous section will degenerate to other forms, such 
that their eigenvalues are still given by the above expressions. 

We choose not to study these degenerate forms since they never appear in the expansion in 
Bell polynomials, however the algebraic dependence of the degenerate A 2n _i on the lower 
degree IOM could be expressed by using Inverse Bell polynomials having as argument the 
tangles D n , of which only a finite number is different from zero. 

We now turn our attention to the A 2n , for this purpose we are interested in the basic 
expressions: 



]Tsm 2 ^ fc (5.5) 



k=l 



for example one can directly check that: 

i 



1!2 2 sin 2 t k = 2L, L > 1 (5.6) 



fc=i 

L 



3!2 4 Y sin 4 t k = 36L, L > 2 (5.7) 

k=l 
L 

5!2 6 Y sin 6 t k = 2400L, L > 3 (5.8) 

k=l 
L 

7!2 8 ^ sin 8 t k = 352800L, L > 4 (5.9) 



k=l 
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which reproduces precisely the coefficients the identity in the A 2n , and in general: 

A 2n = (2n- l)!2 2 "^sin 2rt t fc = 2n[ ^ ~ ^ ) L, L > n (5.10) 



fc=i 



To deal with the eigenvalues of the A 2n -i we will sum up the odd powers of the sine in a 
degeneration independent sum of trigonometric functions, namely: 



£J 3 2 2 - 2 £r V"-*/ sin ((2* -l)i) 



so that our analytic knowledge of the eigenvalues A 2n -i is now completely analytic and ready 
for 1/L expansion, and recovery of the continuum IOM in the leading order behaviour. 
We now observe that by using the following expansion (here B m (z) are the Bernoulli poly- 
nomials) : 

cosqj-l)^) _ 2L , y B 2m (l-l) f (2k-l)n Y^ ( , 

sin (^ - {2 k-l)n + ^ ij r(2m + l)l 2L ) ^ 

and plugging it into the formula for the sum of the odd powers of sines we get: 

2m-l 



7T 2 2 ™~ 3 ^ 2k - 1 \n - k 

j=l k=l v 



r(2m + l)V2L/ 2 2 ™- 2 ^ v 7 V n - 

m=l v ' fc=l v 

(5.13) 

for the contribution of the excitations we have a similar expansion 

y sm— ( t,) = - y ( -ir4- mE ^ (j "' f-) 2m_1 ^r(-i) fc+1 f 2n " X W 

^ KjJ ^ y ' T2m + 1 \2LJ 2 2n ~ 2 ^ y ' \n-k) y 



j£.V m=n x ' k=l 



(5.14) 



however notice that 

(2fc - l) 2m_1 = 0, for m < n (5.15) 



fc=i ^ 



2n - 1 
n — k 



we now define the divergent parts of the eigenvalues A2 n -i as: 



(2n-2)!2 2n - x 7T 2 2n " 3 ^ 2k - 1 V n - k 

y ' k=l 



s 



After introducing the divergent part we can now analyze the finite part of the eigenvalues: 

A 2n _r - Aj^ ~ a m h m -i(V) / M 2 — * J_ " +1 /2n - 1\ 

(2n-2)!2 2 ™- 1 " ^ l J r(2m+l) V2L/ 2 2 ™~ 2 ^ l J \n-kj { 

(5.17) 

being 

Oimhm-l{V) = 2mJ2(j ~ ^) 2m_1 + ^2m.(l " ^) (5.18) 

and 

° n = 2— n!r(3n-I) (5 ' 19) 

where the l2 m -i are the eigenvalues of the Local IOM of CFT, which have been derived in [5]. 
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6 Projectors, Characters and Operator Content 

Let p be the length of the partition V, one has that the projectors: 

project respectively on states with even or odd p. 

It is interesting to mention, that due to the identities: 

B L (A U ...,A L ) _ ((-1) L 1 + R) 

in ~ 2 lb j 

one obtains a TL decomposition of the projectors, which is actually very non trivial, namely, 
for even L: 

n+ = B L (A U A L ) (6 3) 



2L! 
i- 

2L! 



ir = i-%^ (6.4) 



and for odd L: 

B L (A U ...,A L ) 

n = (6 - 5) 

n+ = i + ^ (Ai 2L ," Al) (6.6) 

We then consider the following characters of the projectet space of states: 

V, odd p 



x L J 



S)= <P i&v (6-8) 



■p, even p 

9 



which, by defining the q-Binomial 

/ \ m—x / 



ml l — q l 

can also be expressed in the following Fermionic form: 



X^(g,5) = j:g 2m2+ " 5 (M (6.10) 

m=0 ^ ' a 



xiM) = E « 2ma " m " (2m " 1>5 2m _i < 6 - n » 

m=l ^ 7 9 

alternatively we also have the following Bosonic forms for the finitized characters: 

X L + (q, 5) = ~ ( II^ + 9*"*) + I^ 1 - <? fc - f )) (6-12) 
fc=i fc=i 

1 L L 

X L -(q, S) = 2 ( II^ + ^ f ) - - <f~ty ( 6 - 13 ) 

k=l k=l 

We remark that these finitized characters have been already object of extensive study, and 

are found, in a slightly less symmetric form for example in pp. 

By use the above formulas one can see that the partition functions behave as: 

Z = Tr(n + T M ) ~ Z div (L J M) X L + (q,6) (6.14) 

Z = Tr(n-T M ) ~ Z div (x,L,M) X -(q,S) (6.15) 
q = e**%* (6.16) 



and 

where 
and 



^ - «p (" E ^* 2 " + " E ^^V"" 1 ) < 6 - 17 > 

n=l v y n=l v y 

we now take a look at the expansions of the L = oo characters: 

x +(q,0) = l+q 3 +q 4 +2q 5 +2q G +3q 7 +3q s +±q 9 +5q 10 +6q n +7q 12 +9q 13 +llq u +. . . (6.18) 

X+ (q, 1) = l+g 2 +g 3 +2g 4 +3g 5 +3g 6 +5g 7 +5g 8 +7g 9 +8g 10 + llg n + 12g 12 +16g 13 +. . . (6.19) 
X ™(q,0) = q+q 2 +q 3 +q A +q 5 +2q 6 +2q 7 +3q 8 +4:q 9 +5q 10 +6q u +8q 12 +9q 13 +llq u +. . . (6.20) 
X_ [q, 1) = g 2 +g 2 +q 2 +q 2 +2qz + 2q 2 + 3g 2 +4g 2 +5g 2 +6g 2 +8g 2 +9g 2 +. . . (6.21) 
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we notice that: 

1) = Xo(q) (6.22) 

X-(qA)=Xi(q) (6.23) 

whre xo an d Xi are the respectively the characters of the irreducible h = Q, h — 1/2 
representations of the Virasoro algebra. 

Now, notice that the other 2 characters enumerate the number of partitions of integers into 
distinct parts with either an even (x+) or an °dd (x^) number of parts. We identify such 
characters as the characters of the irreducible (jg)± representations of conformal weight 

X+(q,0) = q-^ X ±(q) (6-24) 

16 

X-M) = g-^XI(?) (6.25) 



16 



so that for even L the corresponding continuum transfer matrix is decomposed into (0) + © 
(jg)- representations whereas for odd L the content is (^) + © (|)_, where the subscripts ± 
are used to specify the parity of the number of fermi modes which we identify with p. 
The states belonging to the (0) + and (|)_ representations, and labelled by a partition V, 
will all be of the form: 

\V) =ij> 2 kp -i ...V_afci-i|0) (6.26) 

1 ' 2 2 1 ' 

with p respectively even or odd. 

The fermion field is expanded in this CclSG clSi 



For the states belonging to the (^ ± ) representation one has instead 

\ = ih i ill . 

16/ 



V)=^. kp ...^ kl 4) (6-28) 



and for the fermion field: 



m = (6.29) 



where the fermi modes satisfy in all cases the following anticommutation relations: 

{^Pn,^m} = <Wm,0 (6.30) 

All the above conformal characters and fermionic decompositions are well known and appear 
in many works, we point out here for example [1] [8] as particularly understandable. 
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7 Conclusions 



In this paper we have obtained a complete decomposition of the transfer matrix for the 
Ising model with periodic boundaries on the periodic TL algebra. The program of building 
a lattice version of the Local IOM of [3] has been completely carried out, and the precise 
relation of the lattice IOM to the continuum IOM is given through a 1/1 expansion whose 
coefficients are proportional to the continuum IOM. Unsurprisingly we find that the operator 
content of the model is classified by the irreducible representations of the A^3,4 CFT. 
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